Constraining models of the large scale Galactic 
magnetic field with WMAP5 polarization data and 
extragalactic rotation measure sources 

Ronnie Jansson^, Glennys R. Farrar^, Andre H. Waelkens^, 
Torsten A. EniSlin^ 

^Center for Cosmology and Particle Physics, and Department of Physics 
New York University, NY, NY 10003, USA 

^Max-Planck-Institute fiir Astrophysik, Karl Schwarzschlld-Str. 1, 85741 
Garching, Germany 

Abstract. We introduce a method to quantify the quality-of-fit between data 
and observables depending on the large scale Galactic magnetic field. We combine 
WMAP5 polarized synchrotron data and rotation measures of extragalactic 
sources in a joint analysis to obtain best fit parameters and confidence levels 
for GMF models common in the literature. None of the existing models provide 
a good fit in both the disk and halo regions, and in many instances best-fit 
parameters are quite different than the original values. We note that probing 
a very large parameter space is necessary to avoid false likelihood maxima. The 
thermal and relativistic electron densities are critical for determining the GMF 
from the observables but they are not well constrained. We show that some 
characteristics of the electron densities can already be constrained using our 
method and with future data it may be possible to carry out a self-consistent 
analysis in which models of the GMF and electron densities are simultaneously 
optimized. 



1. Introduction 

Among the probes of Galactic magnetism, Faraday rotation measures (RMs) and 
synchrotron radiation are among the best-suited for studying the large scale Galactic 
magnetic field. Other probes are either mostly applicable to the nearby part of 
the Galaxy (starlight polarization) or mainly used to study small scale structures 
like dense clouds (the Zeeman effect). While synchrotron radiation has been used 
extensively for studying external galaxies, Faraday rotation measure has been the 
method of choice for probing the magnetic field of our own Galaxy (e.g., Han et al., 
2006; Brown et al., 2007). To most simply probe the large scale Galactic magnetic 
field (GMF) with synchrotron radiation, a full-sky polarization survey is needed, 
at frequencies where Faraday depolarization effects are negligible and where other 
sources of polarized radiation (e.g., dust) are either negligible compared to synchrotron 
radiation or possible to exclude from the data. The 22 GHz K band of the Wilkinson 
Microwave Anisotropy Probe (WMAP) provides such a data set. 

In this paper we constrain models of the large scale Galactic magnetic field 
common in the literature by a joint analysis of the polarization data in the WMAP5 K 
band and rotation measures of extragalactic radio sources. The synergistic advantage 
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of combining these two data sets goes beyond the benefit of achieving a larger 
number of data points; the two obscrvablcs probe mutually orthogonal components 
(perpendicular and parallel to line-of-sight) of the magnetic field. We calculate the 
for a particular choice of GMF model and parameter choice, and use a Markov Chain 
Monte Carlo algorithm to optimize the parameters of each field under consideration. 
Thus we can quantitatively measure the capability of different models to reproduce 
the observed data. 

The content of this paper is structured as follows: we briefly review the theory of 
synchrotron radiation and Faraday rotation and their connection with magnetic fields 
in section 2. Section 3 explains our general approach to constrain models of the GMF. 
Section 4 describes the two data sets used in the analysis. Section 5 details the various 
3D models of the magnetized interstellar medium we have studied. Section 6 give our 
results and a discussion, and is followed by a summary and conclusions. 

2. Observables and input 

2.1. Faraday rotation measures 

As first shown by Faraday, the polarization angle of an electromagnetic wave 
propagating through a magneto-ionized medium rotates in proportion to its 
wavelength squared, 

A6I = RMA2. (1) 
The rotation measure, or RM, in units of rad m~^, is 

where Ue is the total density of ionized electrons, which is dominated by the thermal 
electron density. We follow the convention where a magnetic field pointing towards 
the observer causes a positive rotation measure. 

2.2. Synchrotron radiation 

Relativistic electrons spiralling along magnetic field lines radiate synchrotron radiation 
(Rybicki, G. B. and Lightman, A. P., 1986). For a power law distribution of relativistic 
electrons - commonly called "cosmic ray" electrons, Ucre - characterized by a spectral 
index s, 

ncre{E)dE oc ncre,oE~^dE, (3) 
where Ucrefi is a normalization factor, the synchrotron emissivity is 

l + a 1-., 

j„ CX TlcrefiB^^ U ^ . (4) 

For a regular magnetic field and the above relativistic electron distribution, the emitted 
synchrotron radiation has a large degree of linear polarization, around 75% for a 
power law distribution of clcc;tron with spectral index s = 3. Observationally, the 
percentage of polarization is typically much lower due to depolarizing effects such as 
Faraday depolarization and the presence of turbulent or otherwise irregular magnetic 
fields which depolarize the radiation through line-of-sight averaging. In this paper we 
are only interested in the large scale regular field, and will thus exclusively use the 
polarized components of the observed synchrotron radiation. 
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Figure 1. Overview of the implemented analysis. 



2.3. Input 

As is clear from sections §2.2 and §2.1, the observables we are using are not 
direct measures of the Galactic magnetic field itself, but convolutions of the various 
components of the magnetized interstellar medium (B^rif, and ricre)- The ideal way to 
constrain models of the GMF would thus be to simultaneously and self-consistently 
model the thermal and relativistic electron distributions together with the magnetic 
field. This is beyond the scope of the present work, but should become feasible when 
even larger data sets become available. Instead we use 3D models of Ue and ricre 
present in the literature (detailed in section 5) and do not vary their parameters, 
except in a separate study of the magnetic field scale height's dependence on the scale 
height of the thermal electron distribution in section 6.4. 

3. Method 

Our strategy to model the GMF is straightforward. As illustrated in figure 1, 
starting from 3D models of the distribution of relativistic electrons (ricre), thermal 
electrons (rig) and a model of the large scale GMF parametrized by a set of numbers 
Pa we produce simulated data sets of polarized synchrotron radiation and Faraday 
rotation measures. Generating the simulated data is done using the Hammurabi 
code (Waelkens et al., 2009). We compare the simulated data sets to the observational 
data and calculate x^, as a measure of the quality-of-fit. This value together with po 
is passed to a Markov Chain Monte Carlo (MCMC) sampler that generates a new set 
of GMF parameters, pi, as described in section 3.2. This scheme is then iterated until 
we obtain a Markov Chain that has sufficiently sampled the parameter space. In this 
way we find the best fit parameters and confidence levels for each GMF model, and 
can compare and quantify the capability of different models to reproduce the observed 
data. 
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3.1. Quality- of- fit 

We define in the usual way, 

2 ^2 (data — model) ^ 

i 

with i labcUing the data points (pixels for synchrotron data, point sources for RMs). 
The variance, af, is of utmost importance and obtaining it is a key accomplishment 
of this work. The variance is primarily not experimental or observational uncertainty 
(which will typically be a negligibk^ contribution) , but astrophysical variance stemming 
from random magnetic fields and inhoniogeneities in the magnetized ISM. In the 
case of extragalactic RMs there is also a contribution to the variance from Faraday 
rotation in the intergalactic medium and the source host galaxy. The method we have 
developed to measure from the data is described in §4.1.1 and §4.2.2. This variance 
map contains valuable information about the random magnetic fields. Indeed, when 
studying turbulent and small scale magnetic fields this variance map is itself one of 
the observables that a given model of the random field should be able to reproduce. 



3.2. Parameter estimation 

To minimize for a GMF model and find the best fit parameters we implement a 
Metropolis Markov Chain Monte Carlo algorithm (Metropolis et al., 1953) to explore 
the likelihood of the observed data. We use the approach outlined in, e.g., Verde et al. 
(2003) . For a set of parameters pj we calculate the (unnormalized) likelihood function 
'^{Pj) = e~5^^(^-'^. We then take a step in parameter space to pj+i = Pj + Ap, where 
Ap is set of Gaussianly distributed random numbers with zero mean and standard 
deviation s (the "step length" of the MCMC). If C{pj+i) > C{pj) this new set of 
parameters is accepted, and if £{pj^i) < C{pj) the new step is accepted with a 
probability V = C{pj+i) / C{pj). If a step is not accepted, a new Pj+i is generated, 
and tested. This process is continued until the parameter space has been sufficiently 
sampled. To decide when this has happened we use the Gelman-Rubin convergence 
and mixing statistic, R (Gelman and Rubin, 1992). In the vein of Verde et al. (2003) 
we terminate the Markov Chain when the condition R < 1.03 is satisfied for all 
parameters. To achieve a smooth likelihood surface we run all Markov Chains to at 
least 50k steps. 



3.3. Producing simulated observables from models 

We use the HAMMURABI code developed by Waelkens et al. (2009) to simulate full- 
sky maps of polarized synchrotron emission and rotation measures from 3D models 
of He, ricre ^ud B. The produced polarized synchrotron sky maps (Stokes Q and U 
parameters) are in HEALPixJ format (an equal-area pixelation scheme of the sphere, 
developed by Gorski et al. (2005)) and are calculated taking Faraday depolarization 
into account. The Galactic contribution to the rotation measure (RM) of extragalactic 
point sources can also be calculated. For details, see Waelkens et al. (2009). 



X http://healpix.jpl.nasa.gov 
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3.4- Disk - halo separation 

Significant magnetic fields in galactic halos are observed in external galaxies (Beck, 
2008). However, the magnetic field in our own Milky Way galaxy has proven very 
difficult to study, and a long-standing question in Galactic astrophysics is the nature 
and extent of a magnetic Galactic halo. For instance it is unclear whether the Galactic 
halo field may be completely distinc;t from the magnetic field in the disk, in the sense 
that their topologies are different. In other words, are the disk and halo fields best 
modeled by the same GMF model (with perhaps different best fit parameters) or 
distinctly different GMF models? To address this question we separate the sky into a 
disk part and a halo part, and optimize each GMF model under consideration using 
only data from one part at a time. 

In the case of polarized synchrotron emission it is very difficult to separate the 
contribution of the large scale field from that of smaller, more nearby structures simply 
because the flirx from a particular structure scales as Thus we apply a mask 

on the synchrotron data that covers the disk and clearly defined nearby structures 
such as the Northern Spur (see section 4.1.2). Extragalactic rotation measures are not 
subject to the inverse square effect since the relevent sources are point-like and the 
RM contribution of some volume of the ISM is independent of its distance from us. 

In this paper we will refer to the 'disk' and 'halo' as distinct regions on the sky, 
where the dividing line between the two is given by the set of directions pointing 
towards z = ±2 kpc at Galactocentric radii i? = 20 kpc, i.e., |6| < 4° at / = 0° and 
|6| < 10° at I = 180°. This division is plotted in, e.g., figure 5. 

4. Data 

4-1- Synchrotron radiation 

In the WMAP K-band (22 GHz) the observed polarized radiation is dominated by 
Galactic synchrotron emission. In the WMAP five-year data release the observed 
emission is further analyzed and separated into foreground components caused by 
synchrotron, dust and free-free emission (Gold et al., 2008). Throughout our analysis 
we will use this synchrotron component as our data set for polarized Galactic 
synchrotron emission. The polarized WMAP component was first included, and 
analyzed in the context of the GMF, in the three-year data release (Page et al., 2007); 
subsequent GMF modeling of that data set was done by, e.g., Jansson et al. (2008) 
and Miville-Deschenes et al. (2008). 

The WMAP polarization data is in the form of two HEALPix maps of the Stokes 
Q and U parameters. The resolution is ~1°, with about 50k pixels on the sky. To get a 
more easily interpreted map we can instead plot the polarized intensity P = ^JQ^ -{-V^ 
and with the polarization angle rotated 90° to depict the projected magnetic field 
direction overlaid on the polarized intensity (see figure 2). 

It is clear from figure 2 that emission from the Galactic disk is strong, with 
polarization angles consistent with a magnetic field parallel to the disk. Also present 
is the prominent Northern Spur, a highly coherent structure stretching from the disk 
to high latitudes, and likely the result of the shock between two expanding nearby 
supernova shells (WoUeben, 2007). 

While structures like the Northern Spur can be masked out when studying the 
large scale GMF, smaller features and irregularities in the polarized emission abound 
due to the turbulent nature of the magnetized interstellar medium. Since we are 
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Figure 2. Polarized synchrotron intensity (color) in a MoUweide projection 
with galactic longitude Z = 0° at the center and increasing to the left, overlaid 
with a texture showing the projected magnetic field directions (i.e., the observed 
polarization angle rotated by 90°). Image created using the line integral 
convolution code, ALICE, written by David Larson (private communication). 




Figure 3. Upper panel: sky maps of Stokes Q (left) and U (right), in 8° 
resolution. Lower panel: sky maps of ag and ajj , respectively. 

interested in the large scale regular magnetic field and want to avoid the random field 
in this study we smooth the synchrotron data to 8°and use maps with 768 pixels, as 
shown in figure 3. 
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4.. 1.1. Estimating the variance To calculate for the synchrotron data an estimate 
of the variance is needed that takes into account not just experimental uncertainties 
but fluctuations in the measured emission due to random magnetic fields. The variance 
should de-weight the of pixels in directions of strong random fields or other localized 
sources of polarized emission. 

We estimate the variance in Q (and similarly in U) by 

3 

Here ctq^ is the variance of the ith pixel in the 8°-pixel Q map, whose pixel has the 

value Qi. The quantity qj is a pixel value from the l°-pixcl Q map, where j enumerates 
the iV = 69 l°-pixels that are within a 4 degree radius of the coordinate on the sky 
corresponding to the center of the 8°-pixel Qi. The lower panels of figure 3 show sky 
maps of GQ and au- 

4.1.2. Masking As pointed out in section 3.4; because the synchrotron flux from 
structures such as supernova remnants scale as the inverse square of the distance to 
the object, flux from nearby structures in the disk can exceed emission caused by the 
large scale GMF in the diffuse ISM. However not all parts of the disk may be polluted 
by bright nearby synchrotron emitters. For example, the large "Fan region" , a strongly 
polarized part of the sky around I ~ 120° to 160° and |6| < 15° has been argued to be 
caused by an ordered, large scale (~kpc) magnetic field oriented transverse to our Une- 
of-siglit and with minimal Faraday ck^polarization (WoUcbcn ct al., 2006). Other such 
regions may exist, but it would require more extensive analysis of polarized emission 
at lower frequencies that probe predominantly the nearby part of the Galaxy (due to 
Faraday depolarization obscuring more distant emission) in order to use them in the 
current analysis. 

In this paper we take the simplest approach and mask out the disk and strongly 

polarized local structures such as the Northern Spur. We use the polarization mask 
discussed in Gold et al. (2008), degraded from 4° pixels to the 8° resolution maps we 
use for the Stokes Q and U data. In cases where an 8° pixel covers a region of both 
masked and unmasked 4° pixels, the larger pixel is made part of the mask. This mask 
covers 33.5% of the sky, shown in figure 4, and leaves 511 pixels unmasked. Because 
only a single of these pixels lies inside the 'disk', as defined in section 3.4, we only 
use polarized synchrotron data when fitting data in the 'halo', hence using 510 data 
points for each of the Q and U maps. 

4-2. Extragalactic RM sources 

We use 380 RMs of extragalactic sources from the Canadian Galactic Plane Survey 
(CGPS, Brown et al. (2003)), 148 RMs from the Southern Galactic Plane Survey 
(SGPS, Brown et al. (2007)) and 905 RMs from various other observational efforts 
(Broten ct al, 1988; Clegg et al., 1992; Orcn and Wolfe, 1995; Minter and Spangler, 
1996; Gaensler et al., 2001), for a total of 1433 extragalactic RM sources. The data 
set contains Galactic longitude and latitude (1° and 6°), rotation measure (RM) 
and observational uncertainty {(Tobs)- Typical values of RMs are ± a few hundred 
radians/m^ in the disk and ± a few tens of radians/m^ at large angles from the disk. 
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Figure 4. The mask used for the polarized synchrotron data. The mask covers 
33.5% of the sky. 

The observational uncertainty reported for individual sources tends to be roughly a 
factor of ten smaller than the typical magnitude of the RM in a given direction. 

Some of these 1433 RMs are multiple measurements of the same extended source. 
Including multiple measurements of one source as if they were different sources would 
lead us to underestimate the total variance in the rotation measure of that region. 
To avoid this we round the 1° and 6° of each source to one decimal, and if there are 
multiple RMs with identical coordinates we replace them with a single RM, whose 
variance (c^^j,) is the mean of the o"o5,'s of the individual measurements; the rationale 
for this is discussed in §4.2.2. This procedure removes 103 'sources', leaving 1330. 

4-2.1. Removing ouiliers The measured RM of an extragalactic source is the sum 
of contributions from the diffuse Galactic interstellar medium (the medium we wish 
to study), the intergalactic medium and the RM intrinsic to the source. In addition, 
sightlines that pass through specific structures such as supernova remnants or H II 
regions can receive very large contributions to their measured RM due the high electron 
density and possibly strong magnetic fields pervading these structures (Mitra et al., 
2003). Thus it is expected that the variance in the distribution of RM values is 
significantly larger for sightlines close to the Galactic plane compared to those at 
large angles to the plane, and this is observed (see figure 6). 

In order to remove the sources whose rotation measure is likely dominated by 
either highly localized Galactic or purely non-Galactic contributions, whose variances 
are intrinsically different, we divide the RMs into a 'disk' component and a 'halo' 
component, in identical fashion to the division described in section 3.4. For each of 
these we calculate the mean and the standard deviation of the RM values. We use 
a modified z-scores method (Barnett and Lewis, 1984) and exclude any extragalactic 
sources (EGS) with an RM value three standard deviations away from the mean. 
These two steps are then repeated until no values of RM are more than three standard 
deviations away from the mean of the remaining sample. This process removes 21 EGS 
from the disk and 72 EGS from the halo. 

4-2.2. Estimating the variance As our objective is to study the large scale regular 
field, the variance o%Qg ^ for the zth EGS should ideally account for rotation measure 
contributions of non-Galactic origin and contributions due to small scale random 
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Figure 5. Top: The full sample of 1433 extragalactic RMs, in Galactic 
coordinates (/ = 0° at the center and increasing to the left). Blue squares represent 
negative RM values, corresponding to a magnetic field pointing away from the 
observer. Red circles show positive RM values. The area of the markers scale as 
the magnitude of the rotation measure. The black dotted lines marks our division 
between the disk and halo. Middle: The final sample of 1090 extragalactic RMs 
(559 in the disk, 531 in the halo) used in the quality-of-fit analysis. Bottom: The 
total variance, ct^gs, of the 1090 RMs used in the data fitting. The area of the 
markers scale as <Tegs s-^'l uses the same scale as the RM plots. 
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Figure 6. Histogram of RMs for the disk region (upper panel) and halo region 
(lower panel). Dashed lines demarcate RM values excluded as outliers (see section 
4.2.1). 



magnetic fields. That is, when these contributions to the rotation measure are large, 
'^EGS i should be large, so that the ith EGS has a small weight in the calculation 
of x^- There is more than one way to calculate this variance, and we settle on the 
following scheme to obtain an estimate of the variance: 

, _ ^ (RM.-RM,)^ 

"EGSA — 2^ +"obs,i^ \' ) 

3 

where j labels the sources closest to the ith source, and a\^^^ ^ is the observational 
variance. To get a decent estimate of the variance we require a minimum of A'^ = 6 
sources. For a source in the disk, we use all neighboring sources within dmin — 3° of 
the ith source. If less than 6 sources are found within Qmin we increase the angular 
search radius until 6 sources are included or we reach Bmax = 6° . If we still have fewer 
than 6 sources, we exclude the ith source from our sample. For EGS in the halo, 
where significant changes in the large scale GMF are not expected to occur at as small 
angular scales as in the disk, we instead use drain = 6° and Omax — 12°. 46 sources in 
the disk and 101 in the halo have an insufficient number of neighboring sources to give 
a satisfactory estimate of the variance and are excluded. The remaining sample, which 
will be used in the quality-of-fit calculation, consists of 1090 extragalactic sources and 
are shown in figure 5; together with the corresponding map of aEcs- 

Having explained the procedure for measuring gegSi we now return to the issue of 
how to treat multiple observations of the same source. The correct procedure depends 
on the particular situation. If the source were pointlike and the observations perfectly 
aligned, then the actual variance in the results of the measurements should in principle 
be the same as the mean value of the cr^^Js for the set of observations. However if 
the source is extended and the measurements are not at exactly the same position, 
their variance probes the cfegs from the foreground turbulent ISM as well as from 
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the RM in the source, which is characteristically much larger than the observational 
uncertainty in each individual measurement. Therefore assigning that variance to the 
combined observation in making the fit would incorrectly reduce its pull in the fit 
compared to sources with a single measurement. With the present dataset of RMs, 
there are only a limited number of cases of multiple measurements and it is possible 
to examine them individually. In only one case is the source extended enough that 
the variance in the observations is significantly larger than the mean of the cr^j,^, and 
for that case the (7%Qg value we derive is nearly the same with either (r^j^g assignment, 
vindicating a posteriori our procedure, described in section 4.2. For future large 
datasets this issue will require deeper analysis. If the prescription used here is not 
adequate, an automatic method will need to be developed to handle all cases, which 
is applicable for both extremes. 

5. 3D models of the magnetized interstellar medium 

5.1. Relativistic electron density 

We use a simple exponential model for the spatial distribution of cosmic ray electrons, 

Ccre{r,z) = Ccre,oexp(-r/ft,r)sech^(2;/ft,^), (8) 

with the scale heights hr = 5 kpc and hz = I kpc (following Page et al. (2007)). 
The normalization factor Cere, o is such that for 10 GeV electrons, Ccre{Earth) = 
6.4 X 10~^cm^'', the observed number density of 10 GeV electrons at Earth (Strong 
et al., 2007). The number density for other energies is calculated assuming a power 
law distribution with spectral index p = —3. 

5.2. Thermal electron density 

The best available 3D model of the Galactic distribution of thermal electrons is the 
NE2001 model (Cordes and Lazio, 2002, 2003), based on pulsar dispersion measures, 
measurements of radio-wave scattering, emission measures, and multiple wavelength 
characterizations of Galactic structure. The model has four different components; a 
thin disk, a thick disk, spiral arms, and some local over/underdense regions (e.g., 
supernova remnants). The model is shown in figure 7. The thick disk has a vertical 
scale height of 0.97 kpc. We adopt NE2001 as our baseline model for the thermal 
electron density. We note, however, that the vertical scale height in this model is 
poorly constrained due to uncertainties in pulsar distances. Gaensler et al. (2008) 
estimate the scale height to be 1.83^2 5 ^P*^' ^^'^ (2008) also suggest 

an increase in the scale height by a factor of two in order to avoid an excessively 
large magnetic field in the halo. In section 6.4.1 we investigate the sensitivity of our 
conclusions on the choice of electron density scale height. 

5.3. Large scale magnetic field 

Numerous models of the large scale GMF have been proposed in the literature over 
the last couple of decades. In this paper we test a representative selection of models. 
Most of the considered models were originally proposed to describe either the entire 
GMF or just the disk field. 

All fields considered here are truncated at Galactocentric radius i? = 20 kpc. For 
all models we set the distance from the Sun to the Galactic center to i?„ = 8.5 kpc. 
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Figure 7. The NE2001 thermal electron density model. The figure shows a 
30 X 30 kpc Galactic disk at z = 0. The grayscale is logarithmic. The light 
colored features in the solar neighborhood are a number of underdense regions. 
The small dark speck in one of the ellipsoidal underdensities is the nearby Gum 
nebula and Vela supernova remnant. 

This is the most commonly used distance in the original published forms of the models 
used in this paper. We note that after a few years of confusion on this topic, the most 
reliable estimate of this distance is now 8.4 ± 0.6 kpc (Reid et al., 2009). 



5.3.1. Logarithmic spirals This class of models contains the most common models 
in the literature (see, e.g., Sofue and Fujimoto (1983); Han and Qiao (1994); Stanev 
(1997); Harari et al. (1999); Tinyakov and Tkachev (2002)). 
We define the radial and azimuthal field components as 

Br = B{r,9, z)smp, Be = —B{r,6,z)cosp, (9) 

where p is the pitch angle of the spiral. The function B is defined as 

B{r,0,z) = -Bo{r)cos (^9 + p\n^^ e"!^!/^", (10) 

with (3 — 1/tanp. The pitch angle p is positive if the clockwise tangent to the 
spiral is outside a circle of radius r centered on the Galactic center. At the point 
{rQ,9 = 180°), the field reaches the first maximum in the direction I = 180° outside 
the solar circle. We set the magnetic field amplitude Bo{r) to a constant value, bo, for 
r < Tc, and Bo{r) oc 1/r for r > Tc. With the vertical scale height, zq, the model has 
five parameters: bo, rc, zq, p and ro. 

The above parametrization is usually named bisymmetric spiral (BSS) and 
exhibits field reversals between magnetic arms (see figure 8). By taking the absolute 
value of the cosine in equation (10) we get the disymmetric spiral (DSS). This latter 
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model is often referred to as an axisymmetric spiral in the literature. However, we 

prefer to reserve the term 'axisymmetric' to mean "independent of azimuthal angle". 
A DSS (BSS) field is a logarithmic spiral field that is symmetric (antisymmetric) under 
the transformation 6 ^ 6 + n. 

Another distinction that can be made is a model's symmetry properties under 
z —z, i.e., refiection through the disk plane. We denote a field as symmetric (S) 
with respect to the Galactic plane if B(r, 6, —z) = B(r, 9, z), and antisymmetric (A) 
if B(r, 0,— z) = — B(r, 0,z). This notation agrees with, e.g., Tinyakov and Tkachev 
(2002), Harari et al. (1999), and Kachelriefi et al. (2007); however it confiicts with 
Stanev (1997) and Prouza and Smida (2003). 

We thus distinguish between four different sub-models of the logarithmic spiral 
on the basis of two symmetries; BSSg, BSSa, DSSg and DSS^. 

5.3.2. Sun et al. logarithmic spiral with ring Sun et al. (2008) test a variety of GMF 
models for the Galactic disk using extragalactic RMs. The authors find the model best 
conforming with data to be an axisymmetric field with a number of field reversals inside 
the solar circle. Following equation (9), Sun et al. defines B{r,9,z) = Di{r,z)D2{r), 
with 



where +1 corresponds to a clockwise magnetic field direction when viewed from the 
north pole. 

Since B(r,6, z) does not depend on 9, this is an axisymmetric field in the true 
sense of the word. Sun et al. take the pitch angle to be the average for the spiral arms 
in the NE2001 model, p = —12°; the other parameters are Bq = = 1 /xG, i?c = 5 
kpc, i?o = 10 kpc and zq = 1 kpc. We take these six quantities as parameters to vary 
in our fit. This field model (hereafter named SuuOSd) is plotted in figure 8. 

Sun et al. (2008) also include a separate halo component, detailed in §5.3.3. In 
section 6 we consider the full, composite Sun08 model, as well as the disk and halo 
components separately. 

5. 3. 3. Prouza-Smida halo field The halo component of the SunOS model is a slightly 
modified version of the toroidal halo component proposed by Prouza and Smida (2003). 
The model (PS03 hereafter) consists of two torus-shaped fields, above and below the 
Galactic plane, with opposite field direction (i.e., antisymmetric in z). The magnitude 
of the field is given by 




(11) 



and 



D2{r) 



+1 r > 7.5 kpc 

— 1 6 kpc < r < 7.5 kpc 
+1 5 kpc < r < 6 kpc 

— 1 r < 5 kpc, 



(12) 




(13) 
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Figure 8. Upper left: An example of a bisymmetric spiral (BSS) field model for 
p = —10°, Tc = 10 kpc and ro = 10 kpc. The corresponding disymmetric spiral 
(DSS) looks the same except that there is no field reversal between magnetic spiral 
arms. Upper right: SunOSD, the favored disk model of Sun et al. (2008). Middle 
left: The field proposed by Brown et al. (2007), based on the NE2001 thermal 
electron density model. Middle right: Halo field model proposed by Prouza and 
Smida (2003), and the halo part of the SunOS composite model. Lower left: 
Toroidal disk field proposed by Vallee (2008). Lower right: The magnetic field 
model by the WMAP team (Page et al., 2007). 
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The parameters used in Sun et al. (2008) are — 10 iiG, = 4 kpc, ^;q' = 1.5 kpc, 
— zf^ = 0.2 kpc for 1^1 < and zf = z^i^ = 0.4 kpc otherwise. This model is 
pfotted in figure 8. 

5.3.4- Brown et al. field In Brown ct al. (2007) the authors propose a modified 
fogarithmic spiral model§ (hereafter Brown07) influenced by the structure of the 
NE2001 thermal electron density model with the aim to explain the SGPS rotation 
measure data only (i.e., the fourth quadrant region, 253° < I < 357°). The model 
has zero field strength for Galactocentric radii r < 3 kpc and r > 20 kpc. Between 
3 < r < 5 kpc (the "molecular ring") the field is purely toroidal (i.e., zero pitch 
angle). For r > 5 kpc, eight magnetic spiral regions with pitch angle 11.5° are defined 
with individual field strength bj . The field in the molecular ring and the spiral region 
corresponding to the Scutum-Crux spiral arm is oriented counterclockwise, and the 
remaining regions clockwise. The field strength in region j has a radial dependence 
\Bj \ = bj/r, with the Galactocentric radius, r, in kpc. The vertical extent of the field 
was not considered, as the model was proposed to explain the measured RMs in the 
Galactic disk. The model field is shown in figure 8. 

In our analysis we generalize this model by introducing three free parameters: a, 
which scales the overall magnetic field strengths (bj) used in the original model; rc, 
such that the field strength is constant for r < rc and \Bj\ oc bj/r for r > rc, and an 
exponential vertical scale height zq- 

5. 3. 5. Vallee field In Vallee (2008) and Vallee (2005) the author models the magnetic 
field in the disk as a perfectly toroidal field consisting of concentric rings of width 1 
kpc. The model we consider here, Vallee08, has nine rings between 1 kpc and 10 kpc 
Galactocentric radii, each with a constant magnetic field strength (see Vallee (2008) for 
details). The field is clockwise as seen from the North Galactic pole, except between 
5 kpc and 7 kpc where the field is reversed. In the published model the distance 
between the Galactic center and the Sun is set to 7.6 kpc. For this reason we rescale 
the radial location of the boundaries between the magnetic rings by 8.5/7.6. We note 
that this rescaling still does not allow an entirely fair comparison of the model. The 
model is shown in figure 8. The only two parameters we vary in the fit is a single, 
overall scaling factor for the field strengths and an exponential vertical scale height. 

5.3.6. WMAP field In Page et al. (2007) the authors cite Sofue et al. (1986) and 
Han and Wielebinski (2002) as reason to model the regular GMF with a bisymmetric 
spiral arm pattern. They choose the model 



where ■tp{r) = ipQ + 'tpiln{r/8 kpc) and x(z) ~ xotanh(2;/l kpc), and let r range 
from 3 kpc; to 20 kpc. In Page et al. (2007) the distance from the Sun to the 
center of the Galaxy is taken to be 8 kpc. By fitting model predictions of the 
polarization angle 7 to the WMAP K-band data by using the correlation coefficient 
Tc = cos{2{'yrnodei-7data)), they report a best fit for xo = 25°, ■^o = 27° and V'l = 0.9°. 

§ Details of the model parametrization are not given in Brown et al. (2007), but obtained through 
private communication with the authors, and not reproduced here. 



B{r,(l),z 



'■) = B()[sm')p{r) cosxiz)f + 

cos "tpi^r) cos x(z)0 + sin x(z)z] 



(14) 
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With these parameters rc,rms=0. 76, the rms average over the unmasked sky. This 
loosely wound spiral is plotted in figure 8. 

A few comments about the WMAP model: Firstly, Page et al. (2007) characterize 
their model as a bisymmetric spiral, while in fact it is not (for a fixed radius, |B| has 
the same value at all azimuths, and thus is an axisymmetric field). Secondly, using the 
rms average (as opposed to the arithmetic mean) of Tc leads to an incorrect estimate 
of the best fit parameters, since, e.g., a pixel with rc = —1 (worst possible fit) is 
indistinguishable to = 1 (best possible fit). Finally equation 14 and the above 
quoted best fit parameters differ from the published quantities due to typos || in Page 
et al. (2007). 

Because the authors of Page et al. (2007) only used the polarization angle in their 
analysis no constraints on Bq were found. When performing the model fitting with 
the WMAP field model in this paper we use a four parameter {Bq, xo: V'Oj V'l) model. 
In accordance with Page et al. (2007) no vertical scale height is used. 

5.5.7. Exponential fields To investigate the efficacy of specific model features in 
lowering the total for the optimized model we test a very simple class of GMF 
models, and iteratively add some complexity. We start with a basic axisymmetric 
model defined by an overall field strength Bq, radial and vertical exponential scale 
heights ro and zq, and pitch angle p. As in the case of the logarithmic spirals we also 
test for symmetry/antisymmetry under the transformation z — > —z. We denote these 
models "Es" and "Ea". 

We also test a somewhat more refined model, by introducing a single field reversal. 
We let Bq — » —Bq for r < Tc, and also allow the pitch angle to be different for r < rc- 
Again allowing for different symmetry under reflection in the Galactic plane, we label 
these six-parameter models "E2s" and "E2a". 

The antisymmetry with respect to the Galactic plane of the signs of Extragalactic 
RMs at \h\ > 15° has been used to support the idea that the magnetic field in the 
halo is antisymmetric with respect to reflection through the plane (Han et al., 1997). 
However, the antisymmetry of the signs of RMs seem to hold only for the inner part 
of the Galaxy i.e., |/| < 100°, which may imply that the GMF in the halo is not 
globally antisymmetric. To quantify this, we also investigate a third pair of models: 
"E2no" and "E2so"- These are identical to E2s, a, except that for E2no the only 
field reversal occurs a.t r < rc, z > (North), and a,t r < rc, z < (South) for E2so- 

6. Results cind discussion 

To quantitatively compare the different GMF models introduced in section 5.3 we use 
the per degree of freedom, also known as the reduced x^j 

Xdof = Xtot/{N-P), 

where N is the number of data points and P is the number of free parameters used in 
the minimization. For each GMF model we find Xdof ^'-'^ three separate cases: using 
the complete data set, using only the data in the disk region, and finally using only the 
halo region data. Each model will thus have three separate sets of best fit parameters. 
There is no obvious way to best present the results graphically; in figure 9 we put 
the reduced x^ of the "single" field models (i.e., not the SunOS composite model, but 

II Errata: http://lambda.gsfc.nasa.gov/product/map/dr2/pub_papers/threeyear/polarization/errata.cfm 
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see §6.1.1) for the three cases on three parallel axes, and connect the points for each 
niodcl. For comparison, we also calculate the reduced for the null case {B = 0), 
which serves as a upper bound on the possible range of Xdof^ ^^'^ marked by a red 
bar in figure 9. 

It should be noted that the relative value of the reduced for different models 
is of greater importance than the absolute value for a given model, since the latter 
is sensitive to somewhat arbitrary choices in the calculation of (Tegs ^nd (TQ_^. We 
also recall that in this paper we are not aiming to find the best possible model of the 
GMF, as that would require including random fields and testing a larger set of models. 
Our purpose is to present a new method of constraining GMF models and draw broad 
conclusions regarding GMF models currently used in the literature. 

From figure 9, a few general things are clear. 

• There is no overall tendency of field topologies to simultaneously be a good (or 
bad) fit to both the disk and halo data. This may imply that the Galactic 

magnetic field in the halo has a topologically distinct structure compared to the 
disk, i.e., that the GMF is not a "single" field with the same form but slightly 
different parameters in the halo. 

• The spread of minimized Xdof ^'^^ '^^^^ data is significant; the models that are 

antisymmetric under z —z are all strongly disfavored. 

• The spread of Xdo/ between models when fitting to the halo data alone is much 
less than when fitting to disk data alone. This is to be expected since polarized 
synchrotron data is only used in the halo and is identical for B and — B, making 
for instance a symmetric model indistinguishable from an antisymmetric model. 
The synchrotron data set also has roughly twice the number of data points as the 
RM data set in the halo. 

• For the class of axisymmetric models "E2", the quality of fit to halo data is 
remarkably sensitive to the model's specific symmetry with respect to the = 
plane. The case where only the inner region of the Galactic magnetic field is 
antisymmetric is strongly favored over the completely symmetric or antisymmetric 
model. 

6.1. Fitting models to entire data set 

As noted above, figure 9 suggests that the magnetic structure of the Galactic disk 
is significantly different from that of the halo. To test this hypothesis further we 
investigate the GMF models with the best performance in the disk and the halo. For 
these two models, SuuOSd and E2]sro, we note in table 1 the reduced x^ for the models 
when fitted to the region (disk or halo) it performs best in, as well as the reduced 
X^ for the same model when fitted to the other regions (disk, halo, all sky) when the 
same parameters are used as for the best performing region. It is evident that there 
is a strong tendency for a model to do poorly in regions not used in the optimization. 
Indeed, this effect is so pronounced that when fitted to the disk (halo) data the SuuOSd 
(E2]sro) model predictions for the halo (disk) is even worse than the null case, B = 0. 
We thus conclude that the disk and halo magnetic fields should be studied separately, 
and later combined into a complete model of the Galactic magnetic field; as done in, 
e.g., Prouza and Smida (2003) and Sun et al. (2008). Imposing continuity and flux 
conservation is in general a non-trivial challenge. 
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Figure 9. The reduced fof ^ selection of GMF models. Each model has been 
optimized to fit the data for the full data set, the disk data, and the halo data, 
respectively. In each case a separate set of best fit parameters has been obtained. 
The model names are ordered according to the reduced for the full sky data 
sample. The red bars mark the reduced for the null case, B = 0. 

6.1.1. SunOS composite model The SunOS model warrants some additional comments, 
as it is unique among those considered in this paper in that it has separate components 
for the disk and the halo. The full model is not included in figure 9, since it only 
makes sense to fit a "composite" model to the full sky data set. Instead we note the 
reduced for two interesting cases in table 1. When the full 11 parameter Sun08 
model is optimized to the full data set, the lowest for any model is achieved. Note 
however that both the reduced Xdisfc ^'^d reduced x^^jq are 0.2 larger than the lowest 
X^ achieved when fit by separate components (SuuOSd and E2no)- From figure 9 it is 
clear that the PS03 halo component is a relatively poor model of the magnetic field in 
the halo and is thus forcing the disk component to depart from its preferred parameters 
to improve the fit in the halo. Because of the large number of parameters (11) we 
refrain from calculating the confidence levels for the best fit parameters of the full 
model. Instead we consider the case when the parameters of the disk component are 
fixed to their optimized values (obtained using disk data only) , except for the vertical 
scale height, zq, which is allowed to vary together with the halo field parameters. The 
best fit parameters are summarized in table 1. We note that in combination with the 
halo field the preferred scale height of the disk component is relatively small (~ 0.5 
kpc). The preferred halo field strength is considerable, but as discussed in Sun et al. 
(2008) this is most likely due to an underestimation of the thermal electron scale 
height (see also §6.4.1 for further discussion). 

6.2. Fitting models to disk data 

The best fit field model for the disk data is SuuOSd, the axisymmetric field with 
multiple field reversals presented in Sun et al. (2008). The best fit parameters are 
summarized in table 1 and figure 1 1 . Note that some of these best fit parameters are 
very different from the choices in Sun et al. (2008). Some differences are expected since 
the parameter values were arrived at using slightly different data and a significantly 
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Table 1. Best fit parameters and Icr confidence levels for the best 
fitting models. The 'Region' column refer to what data set was 
used to optimize parameters. The set of reduced has been cal- 

culated from the model predictions of the listed best fit parameters. 
* The full SunOS model with 11 parameters. 
^ The full SunOS model when keeping the parameters of the disk compo- 
nent (except the vertical scale height) fixed to their best fit values, and varying 
the halo parameters only. *Note: the likelihood function of rc for the E2no 
model has a complicated shape, and thus has a less well defined error range (see 
figure 12). 



Model Best fit Parameters Region xjuu X%sk xLi. 



SunOSD 


Bo = l.l±0.1/uG, 
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Figure 10. Left: The best fit SunOSn-typc model for the magnetic field in disk 
of the Milky Way. Right: The best fit GMF field for the halo data, E2no. The 
plot shows the field for z > 0. In the south {z < 0), the field direction in the inner 
region (r < 8.7 kpc) is fiipped. The relative normalization is arbitrary. 



different method of analysis which probes a much larger parameter space. If, e.g., the 
pitch angle is enforced to the originally proposed value of p = —12° the reduced 
optimized for the disk data increases from 1.35 to 1.44, causing the model to lose its 
lead position. 

The runner-up for the best fit disk field is the Brown07 model. It is not surprising 
that these two models perform the best, since in addition to the parameters allowed 
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Figure 11. The la and 2(T confidence levels for the best fit parameters of the 
SunOSo model when fitted to disk data. 

to vary in the fit (three parameters for BrownOT, six for SunOSD) both models have 
several fixed parameters that were originally obtained using extragalactic RM data 
sets similar to the one used in the present analysis. These fixed parameters make a 
fair comparison of the two models difficult. 

6.3. Fitting models to halo data 

The best fit field model for the halo data is 'E2no', by a significant margin. Best fit 
parameters are summarized in table 1 and confidence levels plotted in figure 12, while 
the field itself is plotted in figure 10. 

Because of the limited radial and vertical extent of the electron distributions, 
the contribution to the simulated synchrotron emission and rotation measures from 
regions far from the center of the disk are effectively zero, regardless of the magnetic 
field in that region. We allow the radial and vertical scale height parameters of the 
E2no magnetic field to vary in the MCMC up to 30 kpc and 10 kpc, respectively. 
Figure 12 makes it clear that very large magnetic scale heights are preferred. If the 
assumed rie, ricre are correct, we interpret this to mean that a fairly constant field 
strength is preferred (the spatial extent of which we cannot infer). Alternatively, it 
can be a sign that the scale heights of the electron distributions are underestimated. 
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Figure 12. The la and 2<t confidence levels for the best fit parameters of the 
E2no model when fitted to halo data. 

6.4- Scale heights of electron distributions 

It should be emphasized that synchrotron emission and rotation measures are not 
observables of the Galactic magnetic field alone, but the magnetized interstellar 
medium (B, n^, ricre)- To properly model the GMF a comprehensive and self- 
consistent treatment - where parameters of all three components of the magnetized 
ISM are optimized simultaneously - would be desirable, which falls beyond the scope 
of this paper. 

In this section we briefly investigate how changes in the scale height of electron 
distributions affect the best fit parameters of some GMF models. This is by no means 
an exhaustive study, and is mainly intended to test the versatility of our general 
method of analysis, and to be considered a precursor to future work. 

6.4-1- The vertical scale height of thermal electrons A result common for all GMF 
models under consideration is that constraints on their vertical scale height are 
extremely weak. For the disk, the majority of RMs lie within |6| < 1.5° and no 
useful constraint in the vertical magnetic scale height can be obtained. For the halo, 
constraints on the scale height are also weak. An illustration of this is that a field 
like the WMAP model which lacks a vertical scale height altogether still can achieve 
a fairly good in the halo. 
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The reason the magnetic scale height is poorly constrained is because the 
observables depend on the product of the magnetic field and an electron density (n^ 
for RM and Ucre for synchrotron emission). The models of these electron density 
distributions have their own vertical scale heights and thus significantly reduce the 
RM and synchrotron emission at larger vertical distances from the disk independently 
of whether a significant magnetic field exists there or not. 

Recent work by Gaensler et al. (2008) has shown that the vertical scale height of 
thermal electrons may be significantly underestimated in the NE2001 model, possibly 
by a factor of two. It is thus important to consider the impact of a change in 
the thermal electron scale height on the best fit vertical magnetic scale height. To 
investigate this wc modify the vertical scale height, ho, of the thick disk in the Cordcs 
and Lazio (2002) NE2001 model (keeping fixed the product of the mid-plane density 
and the vertical scale height, which is constrained by pulsar DMs), and fit the E2isro 
model to the RM and polarized synchrotron data in the halo. All magnetic field 
parameters except Bq and Zq are kept fixed to their best fit values with the original 
NE2001 scale height. As seen from the results, which are plotted in figure 13, the 
preferred vertical magnetic scale height and its uncertainty decrease as the electron 
density scale height is increased. For a large electron density scale height, ho > 1-5 
kpc, the preferred vertical magnetic scale height is zq « 1.3 — 2 kpc. Additionally, if 
the vertical scale height of the relativistic electron density tracks ho and this should 
be increased as well, the resulting best fit value of zq would be smaller still. When 
in a future analysis both random and regular magnetic fields are considered together, 
equipartition arguments could be invoked to relate the magnetic field and cosmic ray 
electron density scale heights. 

6.4-2. Radial scale length of relativistic electrons The analysis method introduced in 
this paper can readily be extended from a formal point of view, to include varying 
parameters of the electron densities as well as of the GMF models. We test the 
feasibility of this by letting the galactocentric radial scale length of the relativistic 
electron model vary (but enforcing the overall normalization such that ncre{Earth) 
is at its measured value). Redoing the parameter estimation for the halo model 
E2no with this added parameter, we find the new best fit parameters to be within 
the previously noted confidence levels, and the best fit radial scale length becomes 
Tcre = 6.6 ± 1.5 kpc, which is close to the fixed value rcre = 5 kpc. 

6.5. Discussion 

An obvious but important comment is that a best fit model is not necessarily a 'good', 
or correct, model. And best fit parameters and confidence levels have little relevance 
if the underlying model is a poor depiction of reality. For example, most of the models 
under consideration in this paper have a pitch angle as a free parameter, but the best 
fit pitch angles differ by a factor of a few between models. Hence, if wc lack confidence 
in the veracity of a given model we cannot trust the optimized value of the pitch angle. 

So how do wc achieve confidence in a GMF model? More data, in particular RM 
sources in "empty" regions of figure 5, will strongly reduce the number of GMF models 
that can reproduce the data well. Of special importance are the large gaps of RM 
data in the disk at 0° < ? < 60° and 150° < I < 250°. At present, a few topological^ 
distinct GMF models considered in this paper reproduce the RM data almost equally 
well, while at the same time making significantly different predictions for the rotation 
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Figure 13. The best fit vertical magnetic scale height vs. the vertical scale 
height of the NE2001 thermal electron density of the thick disk when fit to halo 
RMs. The E2isio field model was used, keeping all parameters except Bg and zq 
fixed to their best fit values. The error bars contain 90% of the MCMC sampled 
points. 



measure in the regions of the disk lacking data. For example, in the disk the SunOSn, 
BrownO? and E2s models fit the data best, yet the latter model has only one field 
reversal and the others three; thus even the sign of RM will differ between the models 
for part of the disk. Filling these gaps in the data is thus crucial. In addition, more 
data at high Galactic latitude will enable regions of the sky to be constrained where 
now there are too few RMs to measure a. Pulsar RMs are abundant in the disk, and 
future work will incorporate this data set as well. In addition, polarized synchrotron 
data in the disk can be used if the contamination due to local synchrotron features 
are modeled and subtracted. 

A complimentary approach is to introduce a Xtheory term that incorporates our 

understanding of physics (e.g., V-i? = 0) and information gained through observations 
of magnetic fields in external galaxies similar to the Milky Way (e.g., that pitch angles 
only vary within some specific range). Adding Xtheory to Xtot would be a first step 
toward a fully Bayesian analysis. It would act to disfavor unphysical models that, at 
present, are able to reproduce the current data well. Similarly, using predictions from 
a theory of galactic magnetogenesis, such as dynamo theory, could possibly provide 
powerful constraints on our freedom to build models of the GMF. However, with 
the genesis of galactic magnetism still not fully understood, we believe it prudent 
to study the Galactic magnetic field in the most empirical way possible. Hopefully, 
a firm observational understanding of the magnetic field of the Milky Way will be a 
significant aid in resolving the theoretical question of the formation of such a magnetic 
field. 
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6.5.1. Model-building When building a model of the Galactic magnetic field it is 
essential to be able to quantify the improvement of the model by the inclusion or 
modification of a particular model feature. This is especially true when models become 
more complex than what can be described by a handful of parameters, as we have 
shown in this paper is necessary. The reduced provides such a tool as outlined 
above. In combination with codes to simulate mock data sets and perform parameter 
estimation, the reduced creates a simulation laboratory where one can easily assess 
and compare the pull on of any proposed model feature. 

7. Summary and conclusion 

In this paper we combined the two main probes of the large scale Galactic magnetic 
field - Faraday rotation measures and polarized synchrotron emission - to test the 
capability of a variety of 3D GMF models common in the literature to account for 
the two data sets. We used 1433 extragalactic rotation measures and the polarized 
synchrotron component of the 22 GHz band in WMAP's five year data release. To 
avoid polarized local fc^itures wc applied a mask on the synchrotron data, covering 
the disk and some prominent structures such as the North Polar Spur. 

We developed estimators of the variance in the two data sets due to turbulent 
and other small scale or intrinsic effects. Together with a numerical code to simulate 
mock RM and synchrotron data (the HAMMURABI code, by Waelkens et al. (2009)) 
this allows us to calculate the total for a given GMF model, and then find best fit 
parameters and confidence levels for a given model. This enables us to quantitatively 
compare the validity of different models with each other. 

From applying our method on a selection of GMF models in the literature we 
conclude the following: 

• None of the models under consideration are able to reproduce the data well. 

• The large scale magnetic field in the Galactic disk is fundamentally different than 
the magnetic field in the halo. 

• Antisymmetry (with respect to z ^ —z) for the magnetic field in the disk is 

strongly disfavored. For the halo, antisymmetry of the magnetic field in the inner 
part of the Galaxy is substantially preferred to both the completely symmetric 
or antisymmetric case. 

We find the existence of several false maxima in likelihood space to be a generic 

trait in the models considered. Specifically, for the best-fit disk model in this work, 
the true maximum is far removed from the configuration considered in its original 
form. However, we caution that this optimized configuration should not be trusted 
as its structure is clearly unphysical. We posit a combination of two explanations for 
this fact: firstly, the underlying model may not be sufficiently detailed to explain the 
GMF or may be otherwise incorrect. Secondly, the lack of RM data in large sections of 
the disk may create both local and global likelihood maxima that would not otherwise 
exist. 

From this and other observations we recommend the following regarding future 
GMF model-building: 

• It is necessary to probe a large parameter space to avoid false likelihood maxima. 

• A more complete sky coverage of RM data is necessary to derive a trustworthy 
GMF model. 
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• Optimized parameters are sensitive to the underlying electron distributions, and 

these quantities should be included in the optimization for self-consistency and 

to achieve correct estimates of confidence levels. 

• When possible, combining separate data sets should be done to maximize the 
number of data points in the analysis, to probe orthogonal components of 
the magnetic field, and to ensure the final model's consistency with the most 
observables. 

We also emphasize that by having a well-defined measure of the quality-of-fit, such as 
provided in this work, the efiicacy of any given model feature (i.e., symmetry, extra 
free parameter, etc.) can be quantified, which is essential when building the next 
generation of models of the GMF. 

In the comparatively near future the number of measured extragalactic RMs will 
increase by more than an order of magnitude (see e.g., Taylor (2009); Gaensler et al. 
(2004); Gaensler (2009)) and Planck data will become available, which together should 
allow large parts of the currently valid model and parameter space to be excluded. 

Acknowledgments 

We thank Marijke Haverkorn for fruitful discussions, and Jo-Anne Brown for helpful 
discussions and for compiling the list of extragalactic RM sources. Some of the results 
in this paper have been derived using the HEALPix (Gorski et al., 2005) package. 

References 

Barnett V and Lewis T 1984 Outliers in statistical data Wiley Series in Probability 
and Mathematical Statistics. Applied Probabihty and Statistics, Chichester: Wiley, 
1984, 2nd ed. 

Beck R 2008 in F. A Aharonian, W Hofmann and F Rieger, eds, American Institute 
of Physics Conference Series' Vol. 1085 of American Institute of Physics Conference 
Series pp. 83-96. 

Broten N W, MacLeod J M and Vallee J P 1988 Ap&SS 141, 303-331. 

Brown J C, Haverkorn M, Gaensler B M, Taylor A R, Bizunok N S, McClure-GriSiths 
N M, Dickey J M and Green A J 2007 ApJ 663, 258 266. 

Brown J C, Taylor A R and Jackel B J 2003 ApJS 145, 213-223. 

Clegg A W, Cordes J M, Simonetti J M and Kulkarni S R 1992 ApJ 386, 143-157. 

Cordes J M and Lazio T J W 2002 ArXiv: astro-ph/0207156 . 

Cordes J M and Lazio T J W 2003 ArXiv: astro-ph/0301598 . 

Gaensler B M 2009 in 'lAU Symposium' Vol. 259 of lAU Sym,posium pp. 645 652. 

Gaensler B M, Beck R and Feretti L 2004 New Astronomy Review 48, 1003-1012. 

Gaensler B M, Dickey J M, McClure-Grifiiths N M, Green A J, Wieringa M H and 
Haynes R F 2001 ApJ 549, 959-978. 

Gaensler B M, Madsen G J, Chatterjee S and Mao S A 2008 Publications of the 
Astronomical Society of Australia 25, 184-200. 

Gelman A and Rubin D 1992 Statistical Science 7, 457-511. 



REFERENCES 



26 



Gold B, Bennett C L, Hill R S, Hinshaw G, Odegard N, Page L, Spergel D N, Weiland 

J L, Dunklcy J, Halpcrn M, Jarosik N, Kogut A, Komatsu E, Larson D, Meyer S S, 
Nolta M R, WoUack E and Wright E L 2008 ArXiv e-prints . 

Gorski K M, Hivon E, Banday A J, Wandelt B D, Hansen F K, Reinecke M and 
Bartelmann M 2005 ApJ 622, 759-771. 

Han J L, Manclicstcr R N, Bcrkhuijscn E M and Beck R 1997 A&A 322, 98 102. 

Han J L, Manchester R N, Lyne A G, Qiao G J and van Straten W 2006 ApJ 642, 868- 
881. 

Han J L and Qiao G J 1994 A&A 288, 759-772. 

Han J L and Wielebinski R 2002 Chinese Journal of Astronomy and Astrophysics 
2, 293-324. 

Harari D, MoUerach S and Roulet E 1999 JHEP 08, 022. 

Jansson R, Farrar G R, Waelkens A and Ensslin T A 2008 in R Caballer, J. C D'Olivo, 

G Medina-Tanco and NcUcn, L. and Sanchez, F. A. and Valdcs-Gahcia, J. F., cds, 
'Proceedings of the 30th International Cosmic Ray Conference' Vol. 2 pp. 223 226. 

Kachelriei3 M, Serpico P D and Teshima M 2007 Astropartide Physics 26, 378-386. 

Metropolis N, Rosenbluth A W, Rosenbluth M N, Teller A H and Teller E 1953 
J. Chem. Phys. 21, 1087-1092. 

Mintcr A H and Spanglcr S R 1996 ApJ 458, 194 +. 

Mitra D, Wielebinski R, Kramer M and Jessner A 2003 A&A 398, 993-1005. 

Miville-Deschenes M A, Ysard N, Lavabre A, Ponthieu N, Macias-Perez J F, Aumont 
J and Bernard J P 2008 A&A 490, 1093-1102. 

Oren A L and Wolfe A M 1995 ApJ 445, 624-641. 

Page L, Hinshaw G, Komatsu E, Nolta M R, Spergel D N, Bennett C L, Barnes C, 
Bean R, Dorc O, Dunkley J, Halpern M, Hill R S, Jarosik N, Kogut A, Limon M, 
Meyer S S, Odegard N, Peiris H V, Tucker G S, Verde L, Weiland J L, Wollack E 
and Wright E L 2007 ApJS 170, 335 376. 

Prouza M and Smida R 2003 Astron. Astrophys. 410, 1-10. 

Reid M J, Menten K M, Zheng X W, Brunthaler A, Moscadelli L, Xu Y, Zhang B, 
Sato M, Honma M, Hirota T, Hachisuka K, Choi Y K, Moellenbrock G A and 

Bartkiewicz A 2009 ArXiv e-prints . 

Rybicki, G. B. and Lightman, A. P. 1986 ISBN 0-471-82759-2. Wiley. 
Sofue Y and Fujimoto M 1983 Astrophys. J. 265, 722-729. 

Sofue Y, Fujimoto M and Wielebinski R 1986 Ann. Rev. Astron. Astrophys. 24, 459- 
497. 

Stanev T 1997 Astrophys. J. 479, 290. 

Strong A W, Moskalenko I V and Ptuskin V S 2007 Annual Review of Nuclear and 

Panicle Science 57, 285-327. 
Sun X H, Reich W, Waelkens A and EnfiUn T A 2008 A&A 477, 573-592. 

Taylor A R 2009 ApJ submitted . 

Tinyakov P G and Tkachev I I 2002 Astropart. Phys. 18, 165-172. 
Vallee J P 2005 ApJ 619, 297-305. 
Vallee J P 2008 ApJ 681, 303-310. 



REFERENCES 



27 



Verde L, Peiris H V, Spergel D N, Nolta M R, Bennett C L, Halpern M, Hinshaw 

G, Jarosik N, Kogut A, Limon M, Meyer S S, Page L, Tucker G S, Wollack E and 
Wright E L 2003 ApJS 148, 195-211. 

Waelkens A, Jaffa T, Reinecke M, Kitaura F S and EnBlin T A 2009 A&A 495, 697- 
706. 

WoUeben M 2007 ApJ 664, 349-356. 

Wolleben M, Landecker T L, Reich W and Wielebinski R 2006 A&A 448, 411-424. 



